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Abstract. We study the Hamiltonian equations of motion of a heavy tracer 
particle interacting with a dense weakly interacting Bose-Einstein condensate 
in the classical (mean-field) limit. Solutions describing ballistic subsonic mo- 
tion of the particle through the condensate are constructed. We establish 
asymptotic stability of ballistic subsonic motion. 



1. Introduction 

In this paper we study a physical system consisting of a very heavy quantum 
particle moving through a very dense, very weakly interacting Bose gas. The gas is 
at zero temperature, and it exhibits Bose-Einstein condensation. It is known that 
the speed of sound in an interacting Bose gas is strictly positive. (The dispersion 
law of sound waves, the Goldstone modes of the Bose-Einstein condensate, is linear 
in the magnitude of their wave vector at the origin in wave- vector space; see, e.g., 
[6].) In a limiting regime where the density of the Bose gas is very large, but the 
interactions between atoms in the gas are very weak, the dynamics of the Bose gas is 
harmonic (the Hamiltonian is quadratic in the creation and annihilation operators 
of the Goldstone modes). 

Quantum systems very similar to this one have been studied in |11[ 125] . Under 
appropriate assumptions on the interactions between the particle and the atoms in 
the Bose gas, one-particle states can be constructed that describe inertial motion of 
the dressed quantum particle at a speed below the speed of sound of the condensate. 
Although there are some results on scattering in such systems, asymptotic stability 
or -completeness of this type of motion has not been established yet. In this paper 
we study related problems in the classical or "mean-field" limit. 

One may argue quite convincingly (see, e.g., [13]) that, in the mean- field limit 
where the density of the Bose gas and the mass of the tracer particle become very 
large, the quantum dynamics of the system approaches the classical Hamiltonian 
dynamics of a point particle interacting with a wave medium (in the sense that time 
evolution and quantization commute, up to error terms that vanish in the mean-field 
limit). The corresponding Hamiltonian equations of motion have been exhibited in 
[13] . The purpose of the present paper is to construct a class of solutions of these 
equations describing inertial motion of the particle and to establish asymptotic 
stability of subsonic inertial motion for a suitable class of initial conditions. 
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In the following, Xt £ R 3 denotes the position of the tracer particle and Pt <E R 3 
denotes its momentum at time t. The Bose-Einstein condensate is described, in 
the mean- field limit, by a complex wave field, fit, belonging to a suitably chosen 
Sobolev space of functions on R 3 . This wave field describes low-energy excitations 
of the ground state of the condensate. In a certain limiting regime described in 
|13) . and for a certain choice of boundary conditions of fit at oo, the Hamiltonian 
equations of motion of the system are given by 

(1) MX t = P t , 

P t = -VV(Xt) + [ VW(- - X t ) Re p u 

ifit = -Afi t +/iRe fi t + W(- - X t ). 

Here, V : R 3 —> R is the potential of an external force acting on the tracer particle, 
and W : R 3 —> R is a two-body potential describing interactions between the tracer 
particle and the atoms in the Bose gas. Furthermore, M > denotes the mass of 
the tracer particle, and fi is a parameter related to the product of the density of the 
condensate with the strength of interactions among atoms in the Bose gas, which 
determines its speed of sound, ^JJi. 

One expects the particle motion to exhibit qualitatively different properties ac- 
cording to whether the initial speed of the particle is subsonic or supersonic. For a 
supersonic initial velocity, one expects the particle to decelerate through emission 
of Cerenkov radiation of sound waves into the Bose-Einstein condensate — its speed 
asymptotically approaching the speed of sound in the condensate. For the time 
being, this expectation remains a conjecture that will be analyzed elsewhere. In 
this paper, we construct solutions of the equations of motion describing inertial 
motion of the tracer particle accompanied by a splash in the condensate moving at 
the same velocity. The main result of this paper states that this type of motion 
is asymptotically stable. More precisely, we show that the system QJ of equations 
has "traveling wave solutions" , 

(2) X t = vt, P t = Mv, fit(x) = lv {x - vt), 

describing inertial motion if the speed of the particle, is smaller than or equal to 
the speed of sound in the condensate, and we show that these solutions are stable 
under small perturbations. 

Simplified statements of our main results are summarized in the following the- 
orems; see Section [3] for the full precise statements. The outlines for their proofs 
can be found in Sections 14.11 and 15.11 

Theorem 1.1 (Inertial motion). Consider the system ([T]) ; with W smooth and 
rapidly decreasing, and with V = 0. Fix also a velocity vector v £ R 3 . 

• If \v\ < yjji, then ([TJ has a unique smooth solution of the form (|2]l. with 

| Re 7„(a:)| < (x)^+ , |Im 7 „(x)| < (x)^ 2 ^ . 

• If \v\ = y/JI, then ([1} has a unique solution of the form with 

I Re lv {x)\ < (x)^+, | Im lv (x)\ < (x)(-^+. 
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• If \v\ > y/Ji, then |T]) has no finite- energy traveling wave solutions, except 
for the non-generic case that the Fourier transform of W vanishes on 

{teR 3 \\e + n-(v-i) 2 = 0}. 

Theorem 1.2 (Stability of subsonic inertial motion). Assume that W is spheri- 
cally symmetric, smooth, and rapidly decreasing, and that its Fourier transform is 
nowhere vanishing. Suppose (X t , P t , /3 t ) is a solution to the system Q]) in the case 
V = 0. // the initial conditions satisfy 

• | v | := M" 1 |P | < and 

• j3o — 7„ is smooth and rapidly decreasing, and \\(x) (/3q — 7-u )||l 2 is suffi- 
ciently small (depending on M, fi, \vq\ and W), 

then there exists an asymptotic velocity £ M 3 , with [v^] < ^J~\i, such that 

iM-iPt-Vaol < (l + t)- 3 . 

Moreover, there exists an asymptotic initial position Xq S M 3 such that 

\X t - Voo t-X \ < (l + ty 2 , 

||Re[/3 t -7^0-X o -W)]|U» < + 
||Im[/3 t -7„ 00 (--Xo- ««,*)] IU- £ + 

In other words, the particle velocity converges to some asymptotic subsonic value 
as t /* oo. Moreover, the trajectory of the particle also converges to an asymptotic 
ballistic path, t h- > Xo+Voot, and the field f3 converges to the traveling wave solution 
associated with this ballistic motion. 

In [12) and in [10] (where a more general result going beyond the weak-coupling 
limit is proven) , some of the authors have analyzed particle motion in an ideal Bose 
gas, that is, for \i = 0. In an ideal Bose gas, the speed of sound is zero, hence there 
is only supersonic particle motion. For the models of |12j and |10j . the heuristic 
picture sketched above is vindicated. The particle decelerates to zero speed (= the 
speed of sound), \Pt\ = 0(t~ a ), for some rate a > | possibly depending on initial 
conditions. The condensate wave fit is proven to approach the static solution of 
equation ([]}. 

The existence of solutions describing inertial motion at non-zero velocity shows 
that the static solutions are unstable if \i > 0. However, for an ideal Bose gas, that 
is, fi — 0, the static solutions are stable (|12[ I13j). 

A variety of physics-style results for related problems can be found in J3TJ [51] . In 
[5], the authors prove asymptotic ballistic motion for a forced particle interacting 
with a wave field. However, in their work they assume an independent wave field 
for each point in space. There is extensive mathematical literature on existence, 
nonexistence, stability, and scattering of traveling waves for nonlinear Schrodinger 
and reaction-diffusion equations. For the latter, see for instance [33], or pQ, for 
recent work. For scattering theory in nonlinear Schrodinger equations, see [29] [30] 
[3 HH [TH [20]. For existence and nonexistence results, see [U [TT1 [22]. For results 
on stability of traveling waves, see [HI [3 [9] [26] [3J [16]. See also [23] for a recent 
review of results on traveling waves in nonlinear Schrodinger equations with nonzero 
conditions at infinity. 

1.1. Notations. Here we summarize some notation used throughout the paper. 

• / 9 nieans that there is a constant C depending on k such that / < Cg. 
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• We use the standard notation 

(x) = (l + \x\ 2 )K x = ^-eS n -\ xeR n . 

\x\ 

• Given / : R™ — > C and X S R n , we define the spatial translate f x of / by 

f x (x) :=f(x-X). 

• Given / : R™ — > C, we let F(f) and / denote the Fourier transform of /. 

• We will also use (•, •) to denote the standard L 2 -inner product on R™, 

(/,<?) = / fg, f,g :R"^C. 

In most of this paper, we will only consider three-dimensional space, n = 3. In a 
few instances, we will also use the same conventions in the case n = 1. 
Next, we will use the following notation. 

• Let 5(R") denote the space of smooth rapidly decreasing functions on R™. 

• For any b e R, we let H b (M. n ) denote the standard (inhomogeneous) frac- 
tional Sobolev function space, defined using the norm 

WfWm = / (0 2b \f(0\ 2 ^. 

• Furthermore, when b < 3/2, we can define the corresponding homogeneous 
Sobolev space, H b (W n ), via the norm0 

\\f\\%, = I \e b \fm 2 dt 

In our analysis, we will often encounter functions with different regularities at 
low and high frequencies. These functions will typically belong to the generalized 
Sobolev spaces H b (M. n ), which treat low frequencies like H a and high frequencies 
like H b . Again, we will restrict our attention to spaces with a < 3/2. The spaces 
H b (W l ) can then be described in terms of the norm 

ii/Hh W = / \t\ 2a \f(o\ 2 da+ I \n\ 2b \fm 2 dti. 
j\t\<i j\i\>\ 

The spaces H b and H b represent special cases, with a = and a = b, respectively. 

The following simple observations relate these if ^-spaces, which are somewhat 
technical in nature, to some more intuitive situations. 

• If / e L 1 (R n ) n H b (W l ), then / e H b (W l ) for any a > -n/2. 

• If / e S{W l ), then / e H b {R n ) for any a > ~n/2 and b e R. 

Finally, we are interested in decay properties of fields and particle. In general, 
there are two obstacles preventing decay: 

• Singular behavior of various operators near the Fourier origin — these are 
properties of our system that limit the optimal rates of decay. 

• Lack of spatial decay of our initial conditions and potential. 



Such spaces can still be adequately defined when b > 3/2, but only as tempered distributions 
modulo polynomials. For simplicity, we avoid such technicalities here. 
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We choose our assumptions such that the decay will only be limited by the first 
obstacle mentioned above. For this purpose, we introduce the following definition. 
We say a function / : R n — > C is sufficiently fast decaying iff 

(3) \\(x) N f\\l* = / (x) 2N \f(x)\ 2 dx<^, 

JR» 

for some large enough N. 

Remark 1.3. In this paper, taking N = 6 is sufficient to prove all results. 

2. Basic properties of the equations of motion. 

In this section we discuss some basic properties of the system of equations (fl). 
Moreover, from now on, we will only consider the case V = 0, in which there are 
no external forces acting on the tracer particle. 

2.1. Hamiltonian Structure and Transformations. The system ([1]) is Hamil- 
tonian, with Hamilton functional 



-\V/3\ 2 + m| Re fi\ 2 + W x Re fi 



(4) H(X,P,P) = 7 ±-\P\ 2 + f 
and with the symplectic form 

uj((X,P,(3),(X',P',f3'))=P'.X-P.X' + Im I [3(3' 



In the preceding expressions, we have, for conciseness, adopted the usual practice 
of identifying vector spaces with their tangent and cotangent spaces. 

Moreover, the system has the following invariance and covariance properties: 

(a) Time translation: If (X,P,f3) is a solution of (JTJ) , and if to £ K, then 

X t =X t _ to , P t =Pt-t , fit (a:) = fit-to (*) 

is a solution of (U}. 

(b) Spatial translation: If (X, P, (3) is a solution of ([T]), and if xo el 3 , then 

X t = x + X t , P t = P t , fi t (x) = fi t (x - x ) 

is a solution of ([T]). In other words, ([T]) has translation symmetry. 

(c) Spatial rotation: If (X,P,fi) is a solution of (flj, and if R £ SO(3), i.e., R 
is a spatial rotation, then 

X t = R~ l X t , P t = R- l P t , fi t {x) = fi t {Rx) 

is another solution of (JXJ) , with W replaced by W(x) = W(Rx). In partic- 
ular, if W is spherically symmetric, then ([1]) is covariant under rotations. 

(d) If (X, P, fi) is a solution of (H), and if A > 0, then 

X t =X t , P t = \- 1 P t , fi t {x)=\-ifi t (x) 
is a solution of ([T]), with M and W replaced by 

M = \- 1 M, W(x) = \-*W(x). 



6 D. EGLI, J. FROHLICH, Z. GANG, A. SHAO, I. M. SIGAL 

(e) If (X, P, (3) is a solution of ((TJ, and if A > 0, then 

X t = XX X - H , P t = X^Px-h, h{x) = A-l/3 A - 2t (A-^) 

is a solution of (JT]), with /i and W replaced by 

ji=\- 2 H, W{x)^\-iw{\- 1 x). 

The last two properties can be interpreted as transformations of ([1]) corresponding 
to changes of units. All these properties can be checked by direct computation. 

By rotation and translation transformations ((b) and (c)), we may assume, where 
convenient, that Xq = and Pq = (0, 0, |Po|). Moreover, by appropriate changes of 
units ((d) and (e)), we need only consider the specific case M = 1 and /i = 1. In 
this setting, the equations of motion are 

(5) X t = P t , 

P t = [ \7W Xt Re (3 t dx, 

Jig. 3 

ifit = -Aft + Re /3 t + W Xt . 

From now on, we exclusively analyze equations ([5]) rather than (fTJ. 

Remark 2.1. As the mass M of the tracer particle has been set to 1, the velocity 
v = X and the momentum P coincide numerically. Depending on the appropriate 
physical interpretation, we will use the notations v t and P t , respectively. 

The symmetries of the system ([5]) yield standard conserved quantities, whenever 
such quantities are well-defined: 

• By time translation symmetry, solutions of ([5]) have the energy T~L defined 
in (Ql (with M = /i = 1) as a conserved quantity. 

• By space translation symmetry, solutions of ([5]) have conserved momentum 

(6) V j {X,P,P) = P j + \lva f pdjp, l<j<3. 

• If W is spherically symmetric, then the rotational symmetry property im- 
plies that solutions of ([5]) have conserved angular momentum. 

2.2. Global Well-Posedness. A basic question is whether solutions to the system 
([5]) exist and are unique, for "general" initial data. The following theorem gives an 
affirmative answer to this question. 

Theorem 2.2. Let a € (-5/2, 1/2), b e [1, oo), and W G L^M 3 ) n H b (M. 3 ). If 

PoeM. 3 , Rc /3 e H b a (R 3 ), Im A, G H b a+1 (R 3 ), 

then the system (JS|) has a unique solution 

X t ,P t € C([0,oo);R 3 ), (Re (3 t ,lm (3 t ) e C([0, oo);H b (R 3 ) x H b a+1 (R 3 )), 

corresponding to the initial data Xq, Pq, and /3q. 

The proof of Theorem [221 which employs rather conventional contraction map- 
ping techniques, is deferred until Appendix [A] 
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Remark 2.3. Note that if 

Re ft € H\R 3 ), Im ft, € H\R 3 ), W e L 2 (R 3 ), 

then the Hamiltonian functional % in Q is well-defined. This case is covered by 
Theorem \2. 2\ with parameters a = and b = 1 . 

While Theorem 12.21 shows that P and (3 cannot blow up in finite time, it does 
not exclude the possibility that P and j3 blow up in the limit t f~ oo. However, if 
the Hamiltonian functional is well-defined, then it must be conserved for solutions 
of 12.21 In this case, one is able to not only rule out blowup, but actually establish 
uniform bounds on both quantities. 

Proposition 2.4. Let W e L 2 (R 3 ), and suppose that 

X t ,P t € C([0,oo);]R 3 ), (Re ft,Im ft) € C([0, oo); H\R 3 ) x H 1 ^ 3 )) 

is a solution of ([5])- Then, for any t £ [0, oo), 

(7) \\Pt\ 2 + J Q|V/3 t | 2 + i|Re ft| 2 ) dx < n(X ,P ,(3 a ) + i||T^|| 2 ,. 

Proof. By conservation of the Hamilton functional, 

(8) H(X t ,P t ,p t ) = H(X ,P ,p ) = U . 
It then follows from the definition of T~L that 

\\Pt\ 2 + f (\\^Pt\ 2 + I Re Al 2 ) dx=U - [ W Xt Re /3 t dx 

<Ho + \\W\\ L 2\\Re /3 t \\ L 2. 

The inequality ((T)) follows, since 

\\W\\ L 4 Re ftHi, < i||VK||| 2 + i|| Re /3 t || 2 2 . □ 

2.3. A Change of Variables. From a mathematical perspective, the form ([5]) of 
the equation of motion for j3 is inconvenient because the operator (— A + Re) is not 
complex linear. As a result, Re /3 and Im /3 will have inherently different regularity 
properties in low frequencies. To address this, we make a change of variables, 
introduced in [TO] for the Gross-Pitaevskii equation. In the new variables, the field 
equation gives rise to an evolution that is unitary (on the iJ^-spaces). 
First, let 

U = (-A)*(l- A)"*, 
and define the R-linear operator U r : S(R 3 ) — > S(R 3 ) by 

(9) U r u = J7(Re u) + i(Im u), 

extended appropriately to various Sobolev-type spaces of functions. Note that both 
U and U r are invertible on the appropriate spaces. 
Next, we set 

B t := U- 1 ^ - U-\Ke ft) + i(Im ft). 
Direct computation shows that, by letting 

H= v/-A(-A + l), 
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the equations of motion in terms of Bt take the form 

(10) X t = Pt, Pt= f VW Xt UReB t , \B t = HB t + W Xt . 
Duhamel's principle yields the following expression for B t : 

(11) B t = e-' lHt Bo - i / e-^-^W^ds. 

Jo 

This representation will be advantageous, as it highlights the dispersive properties 
of the system. 

From now on, we use the notation 

(12) H v := H + iv ■ V, veR 3 , 

as this operator will be utilized repeatedly in our analysis. By Fourier transforma- 
tion, we see that the symbol of H v is given by 

(is) h v (o--=\m-vt = \mo-v-i). 

2.4. Symbol Classes. In subsequent decay estimates, we will need to consider 
various differential operators constructed from V, U, H v , etc.; see Section [2~3l 
In particular, we require specific estimates for derivatives of the symbols of these 
operators. Of special importance is how these symbols behave near the Fourier 
origin. As we will do this repeatedly for several operators, it becomes convenient to 
treat them in a systematic manner, depending on their behavior in large and small 
frequencies. As a result, we define the following classes of functions. 

First, let a, b £ R, and suppose / : R™ \ {0} — > C is smooth. We say that / is in 
the class M a iff for every multi-index a, we have 



(i4) |a Q /(OI </,« ,',„, ; ~ ' £e 




These resemble the usual Hormander classes encountered in pseudodiffcrential anal- 
ysis, except, like for the H b -sp&ces, we treat low and high frequencies separately. 
Note the following observations: 

• The map £ i-> |£| a (0 fc ~ a , C G M"\ {0}, is in the class M b a . 

• Given 01,02,61,62 G R, the product of two functions in the classes 
and M. b a 2 , respectively, is a function in the class A^i^i L- 

The above can be verified through straightforward computations. 

Next, we consider similar behavior in only the radial direction. Given a smooth 
function g : (0, 00) — > R, we say that g is in the class O iff for any integer N > 0, 

(15) \g {N) (p)\ < g , N ( P ) b - N , pe(o,oo). 

Unlike for the classes Ai^, here we will need to track the powers of p explicitly, 
hence we only treat the powers of (p) in a unified way. Note the following: 

• The map p ^ (p) b , p £ (0, 00), is of the form O b . 

• Given b\ , 62 G R, the product of two functions in the class O bl and O b2 , 
respectively, is a function in the class O bl+b2 . 

Again, these properties are verified easily. 

Now, for any 9 £ we can define the restriction 

fg : (0,00) ^C, fe(p)=f(p0). 
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In other words, we represent the arguments of / using polar coordinates, and we 
restrict / to a single spherical parameter. We say that / itself is in the class O b iff 
each f e , 9 £ S n ~\ is in the class O b . 

Note the following properties for some basic operators: 

• The symbols £ and |£| for V and |V|, respectively, are in the class M\. 

• Given any b £ R, the operator U b (see Section |2"T3|) has associated symbol 
u b (£) = \£\ b (0~\ which is in the class Mg. 

Of particular interest are the operators H v , v £ R 3 , defined in (|12p . The associated 
symbol, h v , was defined in Q13p . In particular, we have 

h v . e : (0, oo) R, /i„,<?(p) = - « ■ 6), 

for each £ S 2 . Note that when \v\ < 1, we have the uniform lower bound 

ier%(oi = m-v-i\>i-\v\>o, e e r 3 \ {o}. 

Some other properties of h v , listed below, can also be computed directly: 

• If b £ [0, oo), then (h v ) b is in the class M\ b , with the constants in (fl4|) 
depending on \v\ and b. Moreover, if |u| < 1, then (h v )~ b is well-defined on 
R 3 \ {0} and is in the class Mzf ■ 

• If |u| < 1, then for any 6 £ § 2 , 

(16) K,e{p)-\v\{p)- 

• If v £ R 3 and £ § 2 , then h'l e is in the class O . Moreover, the associated 
constants in (|15[) are independent of v. 

3. The Main Results 

In this section, we state the main results of this paper. These pertain to traveling 
wave solutions of (O of the form 

(17) X t = vt, P t = v, p t (z)='y v (x-vt), 

where the function j v is time-independent. These solutions describe uniform motion 
of the particle, along with a splash, 7„, of the Bose gas accompanying the particle. 
The main results of this section, and of this paper, are the following: 

(1) The existence of subsonic and sonic inertial solutions, and the generic 
nonexistence of finite energy supersonic inertial solutions (Theorem 13. ip . 

(2) The stability of subsonic inertial particle motion (Theorem 13. 3p . 

3.1. Inertial Solutions. The first main task is to show that traveling wave solu- 
tions of the form (|17p exist and are unique. Furthermore, we obtain spatial decay 
estimates for these solutions. The precise statements are given below. 

Theorem 3.1. Let v £ R 3 , and let b £ [-1/2, oo). 

• If \v\ < 1, then given any W £ L 1 (R 3 ) fl iJ 6 (R 3 ), there exists a unique 
traveling wave solution l|17p to ([5]), with 

(18) Re 7 , € f| H b+2 (R 3 ), Im 7 „e f| tf a b+2 (R 3 ). 

a>-§ a>-| 

Moreover, if W is sufficiently fast decaying (see ), then for any e > 0, 

(19) \Re 7v (x)\< e (x)- 3+ *, \Im lv (x)\ < e (x)- 2+e . 
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• If \v\ — I, then for any W e L^R") n H b (R 3 ), there exists a unique 
traveling wave solution (|17j) to ([5]), with 

(20) Re 7 „G f) ^ +2 (R 3 ), Im 7 „e f) H b a +2 (R 3 ), 

Moreover, if W is sufficiently fast decaying, then for any e > 0, 

(21) | Re 7t ,(a:)| < e | Im Jv (x)\ < e {x)^ . 

• // \v\ > 1, then for any W G <S(R 3 ) such that W does not vanish on the set 

UeK 3 l i + \t\ 2 -(v0 2 = 0}, 

there do not exist any finite- energy solutions (|17p to ([5]). 
The proof of Theorem 13.11 is given in Section |4] 

Remark 3.2. Theorem \3. 1\ imvlies the following behavior of the profile "f v : 

• If \v\ < 1, both Re j v and Im j v belong to i 2 (R 3 ). 

• // \v\ = 1, only Re j v must belong to i 2 (R 3 ), though Im j v is in H 1 (M 3 ). 

• If \v\ > 1, neither Re j v nor Im j v must belong to i 2 (R 3 ). 

A special case of traveling wave solutions ([T7|) is the case v = 0, that is, the 
stationary solutions of ((5]). In this case, 70 is real- valued and satisfies 

7o = -(A + l)- 1 W. 

This expression implies that 70(2;) decays at infinity like either e - ^' or W , whichever 
decays at the slower rate. 

3.2. Stability of Subsonic Inertial Solutions. Our main aim in this paper is 
to study the stability of the inertial solutions (fTT)) . in the subsonic case (\v\ < 1). 
Roughly speaking, we prove that if the initial speed of the particle is subsonic and 
the initial field is close to a corresponding traveling- wave solution then the velocity 
of the particle converges to an asymptotic (limiting) subsonic velocity and the field 
converges uniformly to a traveling wave corresponding to that limiting velocity. 

A precise statement of our asymptotic stability result is given below. This result 
implies Theorem II .21 

Theorem 3.3. Suppose W G iJ x (R 3 ) is sufficiently fast decaying and spherically 
symmetric, and let 

X = 0, P = v eR 3 , ftefl 1 !! 3 ), 

be initial data for the system (J5]) , such that: 

• I Wo I < 1; *- e - ; the particle is initially subsonic. 

• W is nowhere vanishing. 

• There is some > 0, depending on \vo\ and W, such that 

\\(x) 4 (f3 - j Vo )\\ L 2 < e . 

Then, there is a global solution (X, P — v, /3) to the system ([5]) with the above initial 
data, in the sense of Theorem \2.2\ Furthermore: 
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• There exist an asymptotic velocity Voo G M 3 , with \voo\ < 1 ; and an asymp- 
totic initial position Xq G R 3 such that 

(22) K-Uool <So(l + ty 3 , 

(23) iJft-Woot-Xol <eo(l + 0" 2 - 

• If V(/3 - 7^) G ^(R 3 ), tAen /or any e > 0, 

(24) l|Re(/3 t - 7 ^ + ^ t )|| L »<(l + t)-I ; 

(25) || Im (fi t - 7 £ +, '-*)||l- S (1 + r 

Theorem 13.31 is proved in Section [5l 

Remark 3.4. // one assumes less (but not too little) spatial decay for /3q — j Vo in 
Theorem lS.SH then one can still prove a correspondingly weaker version of Theorem 
\3.Sl in which one obtains a slower rate of convergence from vt to Uqo . 

4. Proof of Theorem 13.11 

In this section, we present the proof of Theorem 13. 11 i.e., we establish existence, 
uniqueness, and decay properties of solutions describing inertial particle motion. 

4.1. Setup and Outline. Much of the strategy for constructing subsonic and sonic 
solutions is the same. We begin the proof by substituting (fl~7| into ([5|), which yields 
the traveling wave equation 

(26) -iv ■ V 7 „ = -A-y v + Re 7 „ + W. 

Making the change of variables G v = U~ 1 ^/ v (see Section [2~3)) and recalling the 
notation (fl"2|) . Eq. ([26]) becomes the complex-linear equation 

(27) H V G V = -W. 

As the transformation from j v to G v is invertible, it suffices to solve this equation 
for G v . 

Now, one must invert H v in order to solve (|27p for G v , and hence j v , in terms 
of W. After Fourier transformation, (|27l) becomes 

(28) K(0G v (0 = -W(Z), 

where h v (£) is the symbol of H v , given in (|13l) . 

Since h v has a (sharp) lower bound |C|" 1 |^-u(C)l > 1 — \ v \: one expects different 
phenomena in the cases |v| < 1, \v\ = 1, and \v\ > 1. To be more specific, consider 
the operator H v = (— A) _1 / 2 iJ„ which has symbol |CI _1 ^t>(0- Observe that: 

• lies in the resolvent set of H v when |u| < 1. 

• lies at the tip of the absolutely continuous spectrum of H v when |w| = 1. 

• is embedded in the continuous spectrum of H v when \v\ > 1. 

It is thus necessary to treat the subsonic, sonic, and supersonic cases separately 
when solving Eq. ([27f for G v . 

If the solution G v exists, then we can determine the regularity of G v and UG V . 
Since Re j v = Re UG V and Im j v — Im G v , this will yield the desired regularities 
for j v in the subsonic and sonic settings. 

With G v , and hence j v , in hand, we proceed to establish another important 
property of inertial solutions: namely that the velocity vt of the particle is constant, 
= v. This is equivalent to showing that P t vanishes identically. To prove this, we 
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turn to the second equation in (|10[) . In particular, we show that the integral on the 
right-hand side necessarily vanishes when Bt(x) is replaced by G v (x — vt). This 
demonstrates that we indeed have inertial solutions. 

The final step is to establish decay estimates for G v and UG V , and hence for 
j v . This analysis focuses primarily on the behavior of G v (£) near the origin in £- 
space, as the strength of the singularity there determines the optimal rate of decay. 
In particular, we expect less decay in the sonic case, as compared to the subsonic 
case, because the symbol h~ x is more singular near the origin (at least along the 
w-direction) in the sonic case than in the subsonic case. Furthermore, since the 
Fourier transform of Im j v is more singular near at the origin than the Fourier 
transform of Re j v , we expect Re j v to have better decay than Im j v . 

Roughly speaking, every power in spatial decay corresponds to one deriv- 

ative in Fourier space, and every such derivative makes the behavior of the Fourier 
transform near the origin worse. This observation results in intrinsic upper bounds 
on the fastest possible rates of decay of -f v . 

Since the best decay rates tend to involve fractional powers of |a;| _1 , we must 
perform a somewhat subtle analysis of those decay rates. The idea is to obtain 
optimal fractional powers by interpolating between adjacent integer powers. For 
this, we derive, in Appendix IB.41 frequency-localized decay estimates, which yield 
precise information for G v restricted to each individual Fourier band |£| ~ 2 fc . When 
restricted to these bands, the aforementioned interpolation process is trivialized. 
We can then piece together these estimates on the individual Fourier bands into 
the desired decay estimates. 

4.2. The Subsonic Case. Here, we assume that |k| < 1. From our assumptions 
on W, we know that W £ H^M 3 ), for every a > -3/2. Since 

Mer^a-M + ieinr 1 ^ 

we can solve for G v £ ff^QR 3 ), as G v = -h^W, since 

A similar estimate shows that 

ll^||^ + 2 < ||G„|| H .+2 < HWHflj < 00. 

This implies the regularity properties (|18p for j v . 

To show that P t vanishes, we make use of ([5]); in particular, 

(29) P t = (Vr*,Re 7^) = Re 

Since W € L^R 3 ) n H b (R 3 ) and b > -1/2, the integral on the right-hand side is 
finite. Moreover, the quantity on the right-hand side within the bracket is purely 
imaginary, hence Pt vanishes identically. 




KI<1, 
KI>1, 



1 1 tmirmW®\ 2d t ■ 



Indeed, on the individual Fourier bands, these interpolation estimates reduce to Holder's 
inequality; see Propositions lB~Tol and fB. 131 



HAMILTONIAN DYNAMICS 



13 



It remains to exhibit the expected decay of j v . Since G v = — H^W, we can 
apply Proposition IB . 10l with f — W and M = H„. Since the symbol h^ 1 of H^ 1 
is in the class M.z\, Eq. (|78|l implies that 

\(xrG v (x)\<\\(x} 4 W\\ L 2, /xe [0,2). 

The above implies the same decay rate for Im j v = Im G v . On the other hand, 
since UG V — —UH~ 1 W 1 and symbol UH~ X is now in the class A4^ 2 , we have that 

\{x)"UG v {x)\ < \\(x) 5 W\\ L 2, M e [0,3). 

This implies the expected decay rate for Re 7„ = Re UG V and proves (fT9|) . 

4.3. The Sonic Case. To understand inertial motion at the speed of sound, i.e., for 
\v\ = 1, is trickier than for |i>| < 1, because the quantity |£| /iu(£) is now singular 
at the origin. However, this singularity only occurs along one "bad" direction, the 
w-direction. Thus, by leveraging the remaining "good" directions, we still obtain 
regular inertial solutions. 

We again solve for G v , as G v = — /z." 1 W', by requiring that the right-hand side 
belongs to an appropriate function space. As | /!„(£) I -1 < |£| -2 in the region |£| > 1, 
we have that 

J\s\>i \ h v{0\ J\t\>i 

For the low-frequency region, we have that 




with a e (-1/2,1/2), 

It remains to show that the integral on the right-hand side is finite. For this 
purpose, we introduce spherical coordinates, (p,0,<fr), on the unit ball, with g 
[0, 7r] measuring the angle from the ^-direction. Then, 

J\i\<\%>{£)\ 2 " Jo Jo Jo |(p)-COS0| 2 

-L So WTW dpd ^ 

where, in the last step, we used that simp < <p, 1 — coscp ~ <p 2 , and (p) — 1 ~ p 2 
on our domain. By reparametrizing p = r cos a and <p — r sin a (or by integrating 
directly), we see that the above integral is finite whenever a > —1/2. 

The above reasoning implies that G v € i^ 2 (R 3 ), and hence UG V £ H%+ 2 (R 3 ), 
for any a > —1/2. From this, it follows that Re ^ v and Im j v satisfy (|20|) . 

To show that P t vanishes, we again take a look at the right-hand side of (|2U1) . 
Like in the subsonic case, it suffices to show that the integral in (|2"9"|) is finite. The 
high-frequency region |£| > 1 is treated in the same fashion as in the subsonic case. 
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For the low frequencies, |£| < 1, we change coordinates, as before: 

< oo. 

To exhibit the expected decay rates, we let S v be defined as in ([751) . i.e., as the 
Fourier multiplier operator with symbol ((£) — v ■ £) _1 . Writing 

G v = -(-A)~ 1/2 S V W, 

we may apply Proposition lB.13l with f — W and M — (— A)^ 1 / 2 . Since the symbol 
of M is in the class Mz\, Eq. flST]) implies that 

\(xrG v (x)\<\\(x) 3 W\\ L2 , M e [0,2/3). 

Moreover, since C/G^ = -(1 - A)~ 1 / 2 S V W, and the symbol (£) of (1 - A)" 1 / 2 is in 
the class Mq 1 , then (I5T1) yields 

\(x)»UG v (x)\ <\\(x) 3 W\\ L 2, M e [0,1). 

The above decay rates for G v and U G v imply (f^Tj) . 

4.4. The Supersonic Case. Finally, if \v\ > 1, it follows from the definition of 
the energy ^ that it suffices to show that Re j v does not belong to L 2 (R 3 ). This 
is easily seen from the following expression: 



i + l£l 2 -(«-£) 2 



Since the denominator 1 + |£| 2 — (v ■ £) 2 vanishes on a two-dimensional surface, 
Re 7„ does not have finite energy under our assumptions on W . 

5. Proof of Theorem 13.31 



In this section, we present the proof of Theorem l3.3l i.e., we establish asymptotic 
stability of subsonic inertial particle motion under suitable hypotheses on the initial 
conditions. 



5.1. Setup and Outline of Proof. By Theorem l3.1[ there exist inertial solutions 
7„, with v g R 3 and \v\ < 1. Next, we define 

G v = U- 1 ^ = -B- X W 

as in Sections 12.31 and 14.11 

By Theorem E2 there is a global solution (X,P,/3) to Eq. (0, with 

X u P t e C([0,oo);R 3 ), (Re ft, Im ft) E C([0, oo); H^R 3 ) x H\R 3 )). 



As in Section 12731 we define 

B t = U- 1 j3 t , BgCQO.oo);^ 1 ^ 3 ))- 
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Recall that we assumed \vq\ < 1. Suppose T > is such that \vt\ < 1 for all 
t S [0, T), that is, the motion of the particle remains subsonic in this time interval. 
It is then reasonable to introduce a field St by setting 

(30) A = te[0,r). 

Note that <5j describes the deviation of the field fit from the inertial solution j Xt . 
As for /3{, we transform <5 t by defining _D t = U~ 1 S t . 

We now cast the equations of motion for the field fit m a more convenient form, 
in terms of D t , as long as the particle motion remains subsonic. Using the equations 
(fTOf for B t and (|27p for G Vt , we derive the following evolution equation for D, 

(31) it> t = HD t -i(d v G Vt ) Xt •v t = HD t + H- 2 VW Xt ■ v t , te[0,T). 
By Duhamel's principle, this equation implies the following integral equation: 

(32) D t = e- im D a - i f e- iH(t - s) H-^W x ° ■ v s ds. 

Jo 

The equation for the particle acceleration, v t , reads 

(33) v t = ( \7W Xt URe D t . 

Jr 3 

At this point, the proof of Theorem 13.31 consists of two main steps: 

(1) The first step is to show that the acceleration it tends to 0, as t — > 00, and 
that the particle moves at a subsonic speed for all times. This will imply, 
in particular, the existence of an asymptotic (subsonic) velocity Voo- 

(2) The second step involves showing that the field f3 t converges to the field of 
the corresponding inertial solution, 'y x ^ +v °° t - 

These steps are carried out in Sections 15.21 and 15.31 respectively. Here, we briefly 
outline the basic ideas underlying our proof. 

For the first step, we will derive an integral equation for v. Since this equation is 
based on the ansatz (|30p and equations (|3"2"|) and (|33|) , one has to impose a bootstrap 
assumption on v that is adapted to its expected decay rate. In particular, this 
assumption justifies (|30|) up to some time T > 0. Starting from (|32[) and (|33[) . we 
derive an integral equation (with memory) 

(34) v t = r{t)- [ M(t-s)v s ds, te[0,T), 

Jo 

where M is a 3 x 3 matrix-valued function. This equation has the general form 
(1 + L)v = r, where L is an integral operator. The goal is to show that 1 + L can 
be inverted, and that r is small in a suitable sense. 

In order to invert 1 + L, we perform a Fourier transform in the time variable. 
Because Lv is defined by a convolution in the time variable of M with v, it becomes 
a product after Fourier transformation to frequency space. Formally, v can be 
expressed in terms of r as follows: 

ii = (1 + M) _1 f. 

By inverting the Fourier transform, we succeed in solving for and controlling v 
in terms of r. We note that all these steps need to be justified rigorously — their 
validity is not apparent, a priori. 

In the proof, we have to accomplish the following tasks: 



1(5 



D. EGLI, J. FROHLICH, Z. GANG, A. SHAO, I. M. SIGAL 



• First, we prove decay estimates for M. Among other things, these estimates 
validate the Fourier transform of (|34|) . as well as the subsequent inverse 
Fourier transform, which yields an expression of v in terms of r. 

• Next, one shows that 1 + M is an invertible matrix. Here, we use the 
spherical symmetry of W to show that M, with respect to an appropriate 
orthonormal frame, is diagonal. Subsequently, we use some special proper- 
ties of M to show that either 1 + M is positive-definite, or Im M ^ 0. 

• Finally, we derive decay estimates for r: 



Combining these ingredients, we will prove the following decay estimate for v: 



Much of the analytical work revolves around proving decay estimates for M and 
r. Here we rely on the following observations: 

• Each matrix element of M and term of r can be expressed in the form 



Jr 3 

where v G M 3 is subsonic, i.e., \v\ < 1. 
• Since \v\ < 1, the phase h v has no stationary points away from the origin. 
In particular, the radial derivative of h v is uniformly bounded from below. 

As a result, one can apply standard integration by parts arguments in the radial 
variable, with each integration by parts resulting in a factor of t . The number of 
integrations by parts one can take depends on the behavior of the integrand near 
the origin in frequency space. This is captured by the exponent I. Consequently, 
the best decay rate one finds for any one of these terms is determined by I. The 
details underlying these decay estimates can be found in Appendix IB. 21 

For some of these terms, one can exploit additional symmetries of the integrand 
to show that there are no singularities at the origin in frequency space. For such 
terms, one can derive an arbitrarily fast polynomial decay in t, because the number 
of integrations by parts one can perform is arbitrarily large. This observation will 
be important in controlling M. For details, see Appendix IB. 2 1 

The analysis sketched above enables us to obtain a strictly better estimate on v 
than the one imposed in the bootstrap assumption. A standard continuity argument 
yields the desired i _4 -decay for v. This decay implies the existence of an asymptotic 
velocity Voo , as well as the rate of convergence of v t to Voo stated in Theorem 13.31 
Moreover, we readily obtain the asymptotic initial position Xq appearing in the 
associated linear trajectory t M> Xq + Voot. 

It remains to carry out the second main step, which consists of proving con- 
vergence of the field. For this purpose, we study the field St using the Duhamel 
formula (|32[) . and we attempt to control each term on the right-hand side in the 
space L°°. The main tools used to carry out this step are dispersive estimates for 
the propagator, e~ ltH (see Section |2~3|) . which are discussed in Appendix IB. 31 

The basic dispersive estimate comes from |19j , where it is shown, roughly speak- 
ing, that if a function / 6 L 1 (M 3 ) is localized in a single Fourier band, then 
||e _ltff f\\t,<*> decays like i~ 3 / 2 . From this basic estimate, one can derive similar 
decay estimates for e~ ttH Mf, where M is a differential operator whose symbol 
belongs to an appropriate class. Just like in the integration by parts estimates of 



K*)<(l + t) 



t > 0. 



\M <(i + <) 



t > 0. 
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Appendix IB.2| the amount of decay in t depends on how singular the symbol of M 
becomes near the origin in Fourier space. This can be analyzed for rather general 
M using interpolation arguments that are trivialized when one considers a single 
Fourier band. We take full advantage of this flexibility when applying such esti- 
mates to several different operators M, as determined by (|32|) and the definitions 
of D and S. 

By piecing together decay estimates localized in Fourier space, we arrive at the 
expected decay estimates, which imply that St — ft — j** decays as desired. Similar 
to Theorem l3.ll since the Fourier transform of Re St is better behaved at the origin 
than the Fourier transform of Im St, we expect better decay for Re St than for 
Im St- 
ill order to complete the proof of Theorem 13. 3[ we must establish decay for 
fit — J^ +Voot , rather than just for St- To accomplish this, we combine the estimates 
alluded to above with a comparison between the traveling waves j Vt and j Voo . This 
last step is relatively straightforward, involving variants of Sobolev estimates for 
the G v 's in terms of W. 

5.2. Behavior of the Acceleration at Large Times. We wish to prove that 
the acceleration of the particle decays in time, provided the initial deviation, 5q, of 
the field from a traveling wave is chosen small enough, as assumed in Theorem 13. 31 
More precisely, we prove that 

(35) \v t \ <e (l + £r 4 , te[0,oo), 



whence the bounds (|22[) and (|23|) . We recall that the constant £0 appears in the 
hypotheses of Theorem 13.31 and must be sufficiently small. 

In order to prove ([35]) . we resort to a bootstrap argument, which proceeds as 
follows. First, we fix a sufficiently small e oc eo, depending on \vq\ and W. For an 
arbitrary but fixed T > 0, we make the bootstrap assumption 

(36) max (1 + t) 4 \v t \ < e. 

0<i<T 

From (|33[) we see that, since ||<5oIIl 2 is assumed to be very small (with respect to 
ll^llij-i); (EH) automatically holds for small T. 

The goal is to show that if (|36[) holds, and if eo ( and hence e) is sufficiently 
small, then we can prove a strictly better estimate 

(37) max (l+t)' 4 \v t \ < \e. 

0<t<T 2 

A standard continuity argument then shows that (and hence (15?]) ') holds for 
all T > 0. This yields (j35j). 

The next step is to derive a relation for it on the interval t € [0,T). Since e is 
small, the bootstrap assumption (|36j) implies that, for any < t < T , 



C-SS) \ Vt \ < \v \+ f \v T \dr <v* < 1, 







with independent of T. Consequently, the ansatz (|30| (and hence (|31|) - (|33|) ) 
remains appropriate for all t £ [0, T). Plugging (|32|) into (|33|) and using the identity 



^In other words, the particle remains "uniformly subsonic". In particular, in various estimates 
that depend on the velocity v* , constants will not blow up. 
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\JW Xt = e- Xt V WW, we find 

«t = r (t) - Rc (^ w ^J q iUe- iH(t - s)+(Xt - Xs) - v H- 2 VW ■ v s ds^ , 

r (t) = Re (VW Xt , Ue~ iHt D ). 

Next, we expand the second term in the equation for v around Vq. Recalling the 
notation (|12[) . we can expand the exponential factor as 

e -iH(t-s) + (X t -X s )-V _ e -iH(t-s) + f^(v- r -v )dT-V+(t~s)v -V 

_ e _i( t _ s ) H „ o +e -i(t-«)iT„ ^ e / < t (« T -oo)<iT-V _ ^ 

From the above considerations we arrive at 

(39) v t = r (t)- [ M(t-s)v s ds, te[0,T), 

Jo 

where r(t) = ro(t) — ri(i) — r2(t), where 

ri(t) = Re (vW, iUe- i(t - s)H <"> ( e />-- Uo ) dr - v - 1)H~ 2 VW ■ v s ds^ , 



ri 



(t) = Re (vW, iUe- i(t - s)H "<> (#~ 2 - H^) ■ v s ds 



Mij(t) 



and where M is the matrix-valued function 

'Re (d Xi W, iUe- itH "°H- 2 d Xj W) t > 0, 
^0 t<0. 

Equation (|39p can be expressed in the form 

(1 + L)v = r, 

where L is the linear operator 

Lh{t) = [ M(t - s)h(s)ds. 
Jo 

The key observation, stated in the following lemma, is that 1 + L is invertiblc. 
Lemma 5.1. There exists a 3 x 3 matrix function K on R satisfying 

• K(t) vanishes for all t < 0. 

• The following decay estimate holds: 

(40) \K(t)\ < w (l + ty 4 . 

• The following identity holds: 

(41) v t =r(t)+[ K(t- s)r(s)ds. 

Jo 

Proof. See Section EXU □ 

In addition, we will also need the following estimates on r. 
Lemma 5.2. For any t 6 [0, T), the functions r^(t) satisfy 

(42) |r„(f)| <w (l + t)- 4 \\(x)%\\ L2 <e (l + t)- 4 , 

(43) Mi)|<we 2 (l + *r 4 , 

(44) \r2(t)\< W e 2 (l + ty 4 . 
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Proof. See Section E2j2 □ 

Applying Lemma 15.21 and recalling the assumptions on So and W in Theorem 
13.31 we obtain 

\r(t)\ <e {l + t)- A + e 2 {l+t)-\ 
Using (|40|) and ([4"I]). we have that 

\M < (£o + e 2 ) / (1 + t - .s)- 4 (l + s)- 4 ds + (e + e 2 )(l + ty 4 
Jo 

<(eo + e 2 )(l+t)-\ 

for all t < T. Choosing e oc Eq sufficiently small, we arrive at the strictly better 
estimate 

\vt\<\{l + t)-\ 

provided £o is chosen sufficiently small. This completes the bootstrap argument 
and hence proves (|3"5j) . 

The bounds ([52} and ([23} follow easily from (|3"5). By (|35]), |u| is integrable on 
[0, oo). Hence there exists an asymptotic velocity, 

Voo = vo + lirn / v s ds. 
Our hypotheses on vq and our bounds on v imply that |t>oo| < 1 and show that 

/oo 
\v s \ds<e (l + t)- 3 , 

which proves (f22|) . 

To prove existence of an asymptotic initial position, Xq, we show that 

Y t =X t - Voot 

has a limit, as t /• oo. For any < t\ < i2, 

r*2 



Y t2 -Y tl = / (v s - Voo)ds. 



This, together with (|2"2"|) . implies that |Yf a — YjJ 0, as ti,t2 co. It follows that 
Yt has a limit Xo, & s i / oo, and 

/oo 
k-«ooMs<e (i + tr 2 , 

which completes the proof of (|23|) . 



5.2.1. Proof of Lemma \5.1\ The main idea is to apply the Fourier transform in the 
t- variable to (|39)l . but some care is required to justify this. We begin by exhibiting 
some rather straightforward properties of M. We note that 

(46) M u -(i)=Rc(C^(0,i|ei(C)- 1 e- it ' l - ( « ) ^ 2 (00^(0), 

p OO 

M«(«)= / e-^R fi <e < W(0,i|eK0" 1 e-^ (e) CK)6-W(0>d* 



Lemma 5.3. T/ie following properties hold. 
• M satisfies the estimate 

(47) |M(£)| < W (1+t)" 6 ,. i€[0,cx>). 
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• M extends to an analytic function in the lower half-plane, II _ . Moreover, 
this extension is continuous on the closure II_ . 

• For any < I < 4 and any z = to — iy G n_ \ {0}, 

I 



(48) 



d_ 

duj 



M(uj - iy) 



<w (i + M)- 2 . 



Proof. The decay estimate (|4"T)) follows from repeated integrations by parts. More 
precisely, we apply Proposition IB . 51 with / = W, v — vq, and N = 6. In particular, 
M is supported on the half-line [0, oo) and belongs to L 2 (M.) n i 1 (R). This implies 
that M can be analytically continued to the lower- half plane II_ , and that it is 
continuous on 1I_ . 

Next, for any z — u) — iy E C \ {0}, with y < 0, Eq. (|4l)|) implies that 

/•OO 

zMaiz) - i / [fte-^Re^W^),!!^^)"^-"^^^- 2 ^^^))^ 



= -iRe(6^(0,i|el(0" 1 ^ 2 (0C J ^(0) 



V 



"v vs;sj 

-izt 



Using integration by parts (see Proposition IB. 5p . one shows that the (i-)integrand 
in the second term on the right-hand side decays nicely, as t / oo. The first term 
on the right-hand side vanishes, since the L 2 -inner product is purely imaginary. 
Thus, we can integrate by parts again, and this yields 

z 2 m 13 {z) = -Rc&w(OAm~ 1 h-?mw(o) 



-izt 



The second term on the right is well defined by Proposition IB. 51 Moreover, by 
Plancherel's theorem, the first term on the right side is bounded by ||W^||?2- 
This proves (gHJ) for I = 0. If I > 0, then 



\z 2 dlM l3 (z)\ = 



OG 



[d 2 t e-^]t i Re , ^mr^^KoimMo)^ 

By (|47p . the integrand is integrable in t. From here, we can integrate by parts 
twice, as in the / = case, and complete the proof of (14"5|) . □ 



Next, we define the auxiliary functions 



< t < T, 



'0 t < or t> 1 , \ nT 

J M(t- s)v s ds t>T. 

From we obtain the equation 

E v = E r — M * E v . 

By definition, E v is compactly supported. Moreover, Lemma l5.3l implies that E r (t) 
is 0(t~ 6 ), as t /• oo. Thus, we can apply the Fourier transform to obtain 

(49) [1 + M(u)]E v (u) = E r {u), lj ER. 

The matrix on the left-hand side turns out to be invertible. 
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Lemma 5.4. Let zeC, with Im z < 0. 

• If He z > 0, t/ier* Im M(z) is strictly negative definite. 

• If Tie z < 0, then Im M(z) is strictly positive definite. 

• If He z = 0. £/ien M(z) is strictly positive definite. 
In particular, 1 + M is invertible everywhere on n_ . 

Proof. Without loss of generality, we may assume that Vq = (0,0, |i>o|). From the 
explicit formula for M in (j46|) and from the assumption that W is spherically 
symmetric, we conclude that M is diagonal. Thus, it suffices to show that each 
diagonal element Ma has the desired sign. 
We define the quantity 

Qiio = mmy'imTiKoior 2 > o, o. 

Note Qi is strictly positive, since W is assumed to be everywhere nonvanishing. 
For any z £ C, with Im z < 0, we have, by (|46]l . that 



M«(«) 



-12/, 



iQi(€){e^" [/l ^ (€)+z] - e"I fc "o(S)-^]}d^di. 



Because the integrand has good decay when t /*■ oo, we can integrate this directly: 



(50) 



Mu{z) 



i 



i 



Next, we write z = cu — iy, with y > 0. If cj = Re z = then (f50"|) yields 



M 4l (z) = 



1 



(MO 
?i(0 



i 



/i« (0 + 



MO 



M£) + y ; 



> 0. 



From the above, we also see that Mu(z) remains strictly positive even when y 
hence the proof is complete in the case where Re z = 0. 
Next, if w 7^ and y > 0, then 



Im Mn(z) = - Im 



= — Im 
2 



Qi(0y 



i 



i 



■ u> 
1 



ui + yi 
1 



\h V0 (0+^\ 2 + y 2 \K {0-^\ 2 + v 2 



0; 



Since ft.„ is everywhere positive, it follows that |/iu (£) + ljJ \ 2 > \h V[) 

w > 0, and \h VQ (£) + w| 2 < |/i Wo (£) — w| 2 if w < 0. It follows that Im Mu(z) is 

negative if w > 0, while Im Ma{z) is positive if w < 0. 

Finally, we examine the limit y \ when w ^ 0. Writing 



Im Mh(lo) = — lim Im 
2 y\0 



MO 



h Vo ((;)-u + yi 



d£ 
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in polar coordinates and recalling the standard formula Im (a =F iO) 1 = ±7r 1 5(a) 
(in one dimension), we can write the integral above as 

im m u (x) = \ I Qi(0*(MO + «)d£ -If Qi(0*(MO - ^ 



2 

where ^4+ is the integral of Qi over the set where h Vo = —u> and A_ is the integral 
of Qi over the set where h Vo — uj. Q 

Suppose first that u> > 0. Then, A + vanishes, because ft,„ is everywhere positive 
on K 3 \ {0}. In contrast, A_ is strictly positive, since Qi > 0, and since h Vo maps 
onto (0, oo). Consequently, in the case lo > 0, we conclude that 

Im M u (uj) = A + - A- < 0. 

On the other hand, if u> < 0, similar reasoning shows that A+ > 0, while A- = 0; 
hence, Im Mu(lj) > 0. □ 

Applying Lemma f5T4l to Eq. (j49]l. we get 

(51) E v {u>) = [1 + Mioj^Erioj) = {[1 + M(w)]- 1 - l}E r {u)) + E r (u). 

Note that (1 + M)~ l tends to 1 at infinity, hence its inverse Fourier transform is 
singular. We therefore prefer to work with (1 + M ) _1 — 1. 

The matrix 1 + M vanishes nowhere on the real line. Moreover, since M tends 
to 0, as w /• oo (see BE])), it follows that [1 + M(iS)\~ x is uniformly bounded from 
above, for all ui G M. □ Thus, the inverse Fourier transform 

p oo />oo 

(52) K{t)= {[l + M(cj)]- 1 -l}e iwt dw = - / M {uj)[1 + M {oj)]- 1 e lult duj 



is well-defined. 

To prove the bound (|40|) in Lemma 5.1, we first estimate 

4 4 

ll(*) 4 ^lk» < E II^MI 1 + ^M]" 1 }!!^ ^ E II^MIU 1 . 

where, in the last step, we have used the boundedness of (1 + M) _1 . Applying 
Holder's inequality and f4"5]) . we find that 

4 



W) 4 K\\ L ~ < E |Kw) a #,M(w)|| x - < oo. 



fe=0 



This proves (|40| . 

Since Lemma 15.41 implies that the integrand in (1521) is analytic on the lower 
half-plane, a contour integration argument shows that, for any y > 0, 



-l 



- 1 )■ e^duj 



K{t) = e yt i 1 + M(w - iy) 

J —oo I 

If i < 0, then, letting y /*• oc and recalling |48|) . one sees that = 0. 
Finally, using (|40|) and ([5"Tj) . we see that 

£?„ = K * E r + E r . 



4 Since ?t„ does not have any critical points on R 3 \ {0}, these integrals are well-defined. 
^The precise upper bound, though, depends on the properties of W. 
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Since Kit) vanishes when t < 0, this yields (|4T|) . 

5.2.2. Proof of Lernma [57B . First, since Xq = 0, we have that 



\X t \ .1 '■' 



(53) '—r< 7 \v T \dT <v* < 1, 

t t Jo 

where u* is as in (|38|) . We use the Parseval identity to write r$ (t) as 

r (t) =Re ^W(0^ it[h{ °-^ ] UDo(0) 

= Re (iCW(0,e- ith «(«[Rr*o(0 + KKO -1 ^"^)) . 
where v — t~ 1 X t . If t < 1, we can apply Holder's inequality to bound 

ro(*) < HWllffiHJoll^. 

Next, for t > 1, we apply Proposition IB.41 with u = 5, and with (l,b) — (1,0) and 
(Z, 6) = (2, —1), respectively. Combining the two cases yields P2"j) . 
For ri, we first expand 

e /.*(o T -«a)«ir-V - 1 = / — e fK v T'- v o)dT'-V^ T 

= f e^^'- Vo)dT '- v [(v T -v )-V]dT 

jnv T/ -v )dT'.w ( r . v \ rfr ^ 



so that, by Parseval's identity, 

/*t y»T 

n(t)=Re / / / ifi(t,s,r, s')ds'd,Tds, 

JO Js Jo 

fi(t, s, r, s') = (VW, [/e-^'-^^oe^t^'-^oJrfr'-v^, . v#- 2 i, s • VW) 

1 /" r 

u = v + / (v T > - v )dr . 

t-s J s 

Note that v is subsonic, since 

|«| < JZ^ + JZ^j |<Mr' < t>„. 
The idea is the same as that for tq. If i < 1, then we can trivially bound 

ri(t) < \\W\\ H i\\W\\ L 2. 
For i > 1, we note that the symbol 

belongs to the class O -3 (see Section El]). Consequently, applying Proposition [H31 
with / = 2, & = —3, and u = u, yields [j 



\fi(t,s,T,s')\ <(l+t-s)- 5 \v s ,\\v s \\\{x) 5 W\\ 2 



^One can apply an argument similar to one in Proposition IB . 5l to conclude even better decay 



for r\ . 
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Similarly, for we can write 

- H- o 2 = J' ^H- 2 dr = -2i jf" H- t 3 (v t ■ V)dr, 
which yields the identity 

r 2 (t) = 2Re f [ (VW , Ue'^^-o H~ 3 (v T ■ V)(<u s • V)W)drds 
Jo Jo 

= 2Re i /* f'{tW(t) , imr^-^^h-'m ■ *r)(£ ■ v.)W(t))dTds. 



Let / 2 (A s, t) denote the quantity within the above double integral. 
When t < 1, we once again have 

T2(t) < \\W\\ 2 L2 . 



In the case t > 1, since 



is in the class O 4 , we can apply Proposition IB. 41 with I = 1, b = —4, and v = vq: 

\f2{t, S ,T)\<{l+t- S )^\\{x)^W\\l 2 . 

By the bootstrap assumption for v, 

lr 2 (*)l < e 2 /* T(l + t - s)- 4 (l + s)- 4 (l + r)- 4 drds < e 2 (l + t)~\ 
Jo Jo 

This completes the proof of (1441) . 

5.3. Asymptotics of the Field, as t tends to oo. In this part of the proof, we 
examine the asymptotics of the field ft. In particular, we wish to prove the decay 
estimates (j2~4")l and (|23|) , and thereby complete the proof of Theorem 13.31 

5.3.1. Preliminary Estimates. We first collect several estimates that are needed 
in the proofs of |24|) and (|25|) . We begin with some L°°-estimates involving the 
operator H^ 1 defined in (p~2|) . 



Lemma 5.5. Let v,v' <E M 3 , with \v\ < 1 and \v'\ < 1, and let f G 5(M 3 ). 

• TTie following uniform bounds hold: 

(54) n^ff-VlU- + ||#-7lU~ <H ll/IU 

WUH^H-'Vfh- + \\H^H^Vf\\ L ~ < KH \\f\\ L *. 

• The following uniform bound holds: 

(55) WUH-'f - UH-'fU- + WH-'f - H-, l f\\ L ~ <,„,,,„,, |v - t/|||/|| i2 . 
Proof. We first note that the symbol, m(£), of UH^ 1 satisfies 

I ( r)\= !£! < f 1 

1 Wl (OII€KO-«-€l ~ \l€|- 2 ICI>i, 

and hence is in L 2 (M 3 ). Thus, we can estimate 

II^-VIU- < \\mf\W < IMMI/IU= < ll/IU- 
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The symbols of the operators if" 1 , UH^H^V, and H~ l H~, x V are also in L 2 (R 3 ), 
and the bounds (|54|) follow. 

Finally, (1551) follows from the above observations and the identity 

(H- 1 - H-^f = i(v - v') ■ H-'H^Vf. □ 

Let G v = U~ lr y v , with »£l 3 and \v\ < 1, as in Section B~T1 For v and v 1 as in 
Corollary [531 (|55|) implies that 

(56) \\U(G V - G V ,)\\ L ~ + \\G V - G v ,\\ L °c < W>M |« - i/|||W|U». 

The next batch of estimates involves decay in time and forms the backbone of 
the proofs of (|24|) and ([25]) . The main technical ingredient is the dispersive estimate 
given in Proposition IB. 81 

Lemma 5.6. Let v <E R 3 , with \v\ < 1. Furthermore, let t <E (0, oo), and let 

feS(R 3 ). 

• The following decay estimates hold: 

(57) || e - ltH /||L~<t^(||V/|| L1 + ||/|| L1 ), 
\\e- itH Uf\\ L °o <*-*||V/|Ux, 

\\e- itH UH- 2 Vf\\ L ~ <t-%\\f\\ L i. 

• The following decay estimates hold for some small e > 0: 

(58) lle-^trVlk- s ^ 1+£ (l|v/|| L1 + ii/iuo, 

l|e- ltH /||L~ S i- 1+£ ||V/|U 1; 
|| e -«*fl--aV/||z- S r^H/llii. 

Proof. For the first and third estimates of (|57|) . we apply the dispersive estimate 
(fTl)) . with the arguments 

M = I, (a\b\a,l)=(o,0,-\,A, 

M = UH- 2 W, (a', b', a, I) = U, -3, ~, Oj . 

For the second part of (|5"Tj) . we apply (J75J to |V|/, with 

M = U\V\-\ (a',b',a,l) = fo,-l,-i, oV 

The proof of (|5"5)) is analogous: to prove the first and third estimates we again 
apply (|74[) . but with the following parameters, 

M = C/ _1 , (a',6',a,0 = (-1,0, -1- e, 1), 

M = #~ 2 V, (a',b',a,l) = (-1, -3, -1 - e, 0). 

The second part of (|58[) follows from an application of ([75]) to |V|/, with 

M=\V\~ 1 , (a',b',a,l) = (-1,-1,-1 -e,0). □ 
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5.3.2. Decay of S. Using the preliminary estimates just established, we are able to 
prove decay of St, as t /~ oo. Since D — C/" 1 ^, we have the identities 

(59) Re S + iUIm S = UD, U' 1 Re S + ilm S = D, 

so that decay estimates on U Dt and Dt will suffice. 

We begin by analyzing UD t . Using the Duhamel formula (|3"1?)) . we see that 

\\UD t \\ L ~ < \\e- itH (Re 5 )\\ L ~ + \\ e - ltH (Ulm 5 )\\l™ 

_ t 

+ [ 2 \\e- iH ^UH- 2 (v s ■ VW)\\ L oods 
Jo 

+ [ \\e- iH ^UH- 2 (v s -VW)\\ L ^ds 



= h + h + h + h- 

To bound I\ and I2, we apply the first and second parts of (|57p and find that 
h < H(||V(Re So)\\ L x + || Re 5 \\ L i) < r*, 
^a^*~ j ||V(lm*o)||£i <^ f - 

Similarly, for J3, we apply the third part of (j5"l 



f— p — 

h<t-%\\W\\ L i [ 2 \v s \ds<t-? / (l + s)- 4 ds<*-i. 
Jo Jo 



Finally, to bound I4, we apply the second part of (1541) , along with the fact that 
e -tH(t-s) jg a b ounc i ec i operator on L 2 (R 3 ), to conclude that 

h < \\w\\ L 2 f \v s \d s < t(i + ty 4 < (1 + 1)- 3 . 

Combining the estimates for through Z4, it follows that 
(60) ||Re5t|| L - < \\UD t \\ L ~ <t~*. 

The proof of decay of D t in t is analogous. Again, by (|3"21) . we have that 
IIAIU- < He-^CZ-^Re <J )||l~ + || e -«"(Im 5 )\\l™ 



[ 2 \\e- iH ^H- 2 (v s -VW)\\ L oods 
Jo 

f \\e- iH ^H- 2 (v s • VTU)|| L ocds 



= Ji + J 2 + J 3 + Ja- 

We bound J\, J2, and J3 using (j55|) : 

Ji < ^ 1+£ (||V(Re 6 )\\ L x + || Re 5„IM < 
J 2 <^ 1+£ ||V(Im5 )|| L1 <i- 1+£ , 
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For J4, we apply ([54]) (like for I4) and obtain 

•h < \v s \ds<(l + t)- 3 . 
Combining these bounds yields 

(61) ||im<y LO = <||A|U- <^ 1+£ - 

Bounds (ISU)) and (joTj) prove the desired decay for 6. 



5.3.3. Completion of the Proof of Theorem \3.3\ In order to complete the proofs of 
(IM1) and (j2"5)l , we first control the intermediate quantity (3 t — 7^* . Recalling 
and applying (|56p . we obtain that 

(62) || Re [fit - T^)||i- < || Re ^Hi- + || Re ?7(G„ t - G„JIU- 

<t-^ + \v QC -v t \\\W\\ L 2 

/OO 
|v s |ds 

By similar reasoning, using (1561) and (|61[) . we also have that 

(63) ||Im(A-7S)IU~ <*~ 1+£ - 

Finally, in order to establish (IM1) . we use the triangle inequality 

II ^ {fi t - tS^-^IU- < II Re (A " T&)||l- + || Re ( 7 ^ +,w - 7^)IU- 

The first term on the right-hand side is controlled by using (j6"2")l . To bound the 
remaining term, we apply the first inequality in ((BT)) to VW and then use (|23p: 

l|Re(75 +1,oot -7S)llL~ < H^GS^-G^JIU- 

< |X + Woo i-X t |||C/VG Woo || L oo 

< yi/fl^vwiiL-cn-t)- 3 

<||VVF|| L2 (l+t)- 2 . 

These estimates result in the bound (|24p . The remaining bound (|2"5j) for the imag- 
inary part of the field is proven in a very similar manner. 
This completes our proof of Theorem 13.31 



Appendix A. Global Well-Posedness 

In this section, we establish the well-posedness properties of the system ([5]). 
More specifically, here we will prove Theorem 12.21 The proof is an application of 
standard methods, in particular the contraction mapping theorem. The process 
here is particularly straightforward, as this system is close to linear. 
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A.l. Local Well-Posedness. The first step in the proof of Theorem 12.21 is the 
following general local well-posedness result. 

Proposition A.l. Fix a G (—00, 1/2) and b G [l,oo), let Po G R 3 , and suppose 
both ipo € H b (R 3 ), 7To G P^ +1 (R 3 ) are real-valued. Then, there exists T > 0. 
depending on \\W\\ H b ^, \\(po\\ H b, and \\tto\\h'> y > such that the system (J5]) has a 
unique solution 

A,P e C([0,T];R 3 ), (Re /3,1m (3) G C([0,T}; H b a (M 3 ) x H b a+1 (R 3 )), 
with corresponding initial data X$ = 0, Pq, and (ip o,7To) . 
Proof. For convenience, we let 

n = H^ll^ +1 ' R = W^Whi + lkoll^ +1 , 

and we make the usual change of variables, B = U^ 1 ^, as in Section l2~3l Integrating 
the second equation in (fT0|) and applying (fTT|) . we obtain the fixed point equation 
<1>(P) = P, where the map from C([0,T];R 3 ) into itself, is given by 

(64) MQ)]t = Po +Re / I (VW) Ys e- iHs UB dxds 

Jo Jr 3 

-Rei f f (VW) Y ° f e - iH{s - T) UW Y -dTdxds, 
Jo Jm 3 Jo 

with Y G C([0,T];R 3 ) defined 

Y t = X Q + [ Q s ds = [ Q s ds. 



It suffices now to show $ is a contraction on C([0,T];R 3 ) for small enough T. 
Given Q, Q' G C([0, T]; R 3 ), we have 

(65) [$(Q) - $(Q')]t = Re / / h{s) ■ e- iHs UB dxds 

Jo Jr 3 

-Rei / / / / 2 (s,r) ■ e-' lH{s - T) UWdxdTds, 
Jo Jo Jr 3 

/ x ( s ) = (W) Ys - (VW) r » 

7 2 (s,t) = (W) y -- y - - (VW)^'- r -, 

where V is constructed from Q' in the same manner as Y from Q. To bound (|65l) . 
we begin by writing 

||3(Q) - $(Q')IU~ < T sup ll/i^l^-.He-^L/Poll^ 
«e[o,T] 

+ T 2 sup \\I 2 ( S ,r)\\ H - b \\e- iH ^UW\\ Hb 

s,t£[0,T] 

<T\\B \\ Hb sup ||7i(s)|| H - 6 
a+1 se[o,r] 



+ T 2 ||W|| H * sup \\I 2 (s,t)\\ h -j 



1+1 s,t£[0,T] 



<TP sup \\h{s)\\ H - b +T 2r R sup ||I 2 (s,r)|| ff - 

se[0,T] a s,t£[0,T] 
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As a result, our goal is to control Ii(s) and his, t). By the mean value theorem 
and the translation-invariance of the i^-norms, we obtain 

sup < \\Y -Y'\\ L ~\\V 2 W\\ H -» 

se[o,T] 

<T\\W\\ H ->+2\\Q - Q'\\ La o 
<TK\\Q~Q'\\ L ~. 

In the last step, we used the assumptions b > 1 and a < 1/2 to bound the HZa+2' 
norm by the iJ^ +1 -norm. Similarly, we can bound I2{s,t) as follows: 

sup \\I 2 (s,T)\\ H - b < \\Y-Y'\\ L ~\\V 2 W\\ H - b <TK\\Q-Q'\\ L ~. 

s,re[0,T] 

Finally, combining all the above, we obtain 

||$(Q) - $(Q0IU- £ K{T 2 R + T 3 7l)\\Q - Q'|| £ «. 

Therefore, if T is sufficiently small, depending on i? and 71, then $ is a contraction 
on C([0,T];R 3 ). We set P to be the unique fixed point of $. 
With P in hand, we can now define X by 

X t = I P s ds, 



and then B using the Duhamel formula (fTTj) . In particular, since VF has finite 
^ +1 -norm, then (TTT]) implies I? € C*([0, T]; i?^ +1 (R 3 )). It is clear that X, P, and 
B form the unique solution to (ITU)) in the given spaces, with initial data X n = 0, 
Pq, and Bo, respectively. Finally, defining f3 t = U r B t for each t g [0, T], then 

Re G C([0,T];ff Q b (R 3 )), Im = C([0, T}; H b a+1 (R 3 )). 

Moreover, X, P, and j3 must form the unique solution to <j5j) in the desired spaces, 
with initial data Xq =0, Pq, and (</3o,7ro)j respectively. □ 



A. 2. Global Well-Posedness. We can now complete the proof of Theorem 12.21 
First of all, since W G L 1 (M 3 ) n H b (R 3 ) by the assumptions of Theorem [2~2| then 
W G H b +1 (R n ) for all a G (-5/2, 1/2). Suppose we have a local solution of ©, 

X,Pe C([0,T);R 3 ), (Re/3,Im/3) G C([0, T); iJ Q fc (M 3 ) x i/ a b +1 (M 3 )), 

for a fixed a in this range, with T < oo. Since the time of existence from Proposition 
I A. II depends only on the size of the initial data for j3, then it suffices to show 

(66) sup \\B t \\ H t (R3) < C, 

0<t<T a+1 

where B t = U, r T 1 f3t- If holds, then one can apply Proposition I A. 1 1 with initial 
data (Xt-c, Pr-e, Pt-c), with sufficiently small e > 0, in order to push the above 
solution beyond time T . This implies that there cannot be a finite maximal time 
of existence and hence complete the proof. 
To establish (|66|) . we can simply use ([TTT) : 



\B t \\ Hb <\\e- itH B \\ Hb +T sup We-^-'^W 
a+ ° + «e[o,T) 

<p |Uj +i +T||W|| H , +i . 



A" 3 



This concludes the proof of Theorem 12.2 
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Appendix B. Decay Estimates 

We collect in this appendix the decay estimates, in both time and space, that 
are used in the paper. First, we give a brief overview of Littlewood-Paley theory, 
which will be used in the decay estimates of Sections IB. 31 and IB.4I 

B.l. Littlewood-Paley Theory. Let ip e <S(R 3 ) be compactly supported in the 
annulus 1/2 < |£| < 2. For each k G Z, we let tpk G S(R 3 ) be its rescaling 

In addition, we can choose ip such that 

5>*(0 = L, £eR 3 \{0}. 

feez 

We can also assume ip is spherically symmetric. 

Remark B.l. Note that by scaling properties, ipk is in the class A^g. 

Recall that the usual (homogeneous) Littlewood-Paley operators arc defined as 

P fe : S(R 3 ) -> S(R 3 ), P^u = t/j k u, k e Z. 

In other words, P^u is a "projection" of u to the frequency band |£| ~ 2 k . For 
convenience, we sometimes adopt the notation 

P~k = Pk-l + Pk + Pk+l, fcGZ. 

Due to Fourier support considerations, we have Pk = PkP~k- 

Next, we recall the Bernstein, or finite band, inequalities, cf. 0132], which use 
the operators Pk to convert derivatives to constants, and vice versa. 

Proposition B.2. Let p e [l,oo], and let f e S(R 3 ). 

• IfbeR, then 

(67) ll^|V| b /l|L P <2 bfe ||/|| iP . 

• In addition, 

(68) \\PkVf\\ LP <2i/|| LP; \\P k f\\ L , <2- fc ||V/|| iP . 

In this paper, we require more than the usual finite band estimates found in 
(|67p and {55]). We need extensions of these estimates to other Fourier multiplier 
operators. Examples of such operators include U and powers of H v , where set 3 . 
The generalized estimate we need is given in the subsequent proposition. 

Proposition B.3. Let a, b G R, let p G [l,oo], and let M be a Fourier multiplier 
operator, with symbol m : R 3 \ {0} C in the class A4^. Then, for any f £ S( 



(69) \\PkMf\\ LP 



< 



2 bfc ||/||Lp k>0, 
2 afe ||/||L P £;<0. 



Proof. The proof is a straightforward adaptation of standard Fourier multiplier 
estimates; see, e.g., Lemma 2.2]. □ 
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B.2. Subsonic Decay Estimates. First, we discuss a simpler class of stationary 
phase estimates, in which the radial derivative of the phase never vanishes. In this 
case, we can always integrate by parts in the radial direction. 

In particular, this situation holds when the phase is the symbol h v of H v and 
when \ v\ < 1; see (|16[) . We encounter precisely this situation in the proof of Theorem 
13.31 in particular for controlling the acceleration of the tracer particle. 

Proposition B.4. Assume the following: 

• Let I eZ and b <G R, with I > and b < 0. 

• Let v e M 3 . with \v\ < 1. 

• Suppose m : M. 3 \ {0} C is in the class O b . 

• Let f,g € S(R 3 ). 

Then, for any t € (0, oo), we have the estimate 



(70) 



L 2 ■ 



Proof. Throughout, we write O r to denote any function in the class O r . Expressing 
(f70f in polar coordinates and letting F = fg, it suffices to derive 



s 2 Jo 



- ith «.»(.P)p l + 2 O b (p)Fg{p)dpde 



<t- l - 3 \\{x) l+s f\\ L 4{x) l+3 g\\ 



L ' 2 ■ 



Let I denote the left-hand side of the above inequality. 
Recalling (|16p . we bound 

_d_ e —ith V: e(p) 



L = r 



< f 



S 2 JO 



P l+2 o\ P y 



-F e {p)dpd6 



I r p l+2 O b - 1 (p)^e- lth ^^F e (p)dpdd 
Js 2 Jo d P 



Integrating the above by parts yields 



/ < r 



s 2 Jo 



e -ith v , e (P) b-l {p) [p l+l Ff) {p) + p l+2p, {p)]dpd0 



Moreover, we can iterate this process until all positive powers of p are exhausted: 



L < t- 1 - 3 



s 2 Jo 



\O b - l - 3 {p)\\F {p)\dpde + t- l - 3 \F B {Q)\ 



+ t 



-1-3 



s 2 Jo 



1+2 



\o b - l -*(p)\Y,p k \Ft +1) (p)\dpde. 



k=0 



Note that in the final integration by parts, since no powers of p remain, we pick up 
an extra boundary term at the origin. 

Noting our assumptions for I and 6, we have 

I Z *~'" 3 (II*1U~ + llVFHi- + ||V 2 F|| L oo) 

i 



+ t 
<h+h 



-1-3 



s 2 Jo 



O b - l -\p)p^P k \H k+ "\p)\dpd9 



k=0 
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For I\ , we can trivially bound 

h<t~ 1 - 3 E II^/IU-II^IU-^IK^VII^IK^II^. 

\a\ + \P\<2 

For I2, we convert back to Cartesian coordinates to obtain 

h < t-^ 3 E \\d a fh4d l3 g\\Li < t- l - 3 \\(x) l + s f\\ L 4(x) l + 3 g\\ L ,. 

\a\+\0\<l+3 

Combining the above completes the proof. □ 

Proposition B.5. Let v g R 3 , with \v\ < 1, and let f E S(M. 3 ), with f an even 
function. Then, for any N > 3 and t € (0, 00 ), we have 

(71) 

Proof. For 1 < i, j < 3, let 

iij-Im (dJ,Ue- itH »H- 2 d 3 g). 

By Holder's inequality, we can trivially estimate Ii.j < ||/||j,2- 
Writing Zy in polar coordinates, we have 

"°° p*\f{p6)\ 2 e- H P^- B ^ 



lm(Vf,Ue- itH "H- 2 Vf)< HJ , N (1 + t)~ N \\ (x) N f[ 2 



L- ■ 



Iij — Im 



S 2 



H u 



{p){{ P )-e-vf 



-dpd6»=:im / e l e j i e de. 



Fixing 9 £ S 2 , we see there are smooth even functions 



Ug = 



Ug, Wg : K — > K. 

such that 

/>oo 

h = / p 3 ug{p)e' itpw ^Ap 



\f{pQ)t 



(p)((p)-v9)i> 



u>e(p) = {p)-v- 



_ 1 
~~ 2 
1 

~~ 2 

where we applied a change of variables p — > —p in the last step. Thus, 

1 r°° 

' 1 p 3 ug(p)e- u ^-^dp, 
and recalling the explicit form of ug, we have 



p 3 ug(p)e' itpW(, ^dp - J p 3 ug(p)e itpwe{p) dp 
p 3 ug(p)e- itpw »^dp+ / p 3 ue{p)e- itpw o {p) dp 



Oj / p 3 0- 3 {p)\f{ P 6)\ 2 e- itp ^- v ^dpd9. 



From here, the proof proceeds like for (|70|) , except we no longer have a boundary 
at p = 0. Because of (fT6|) . we can integrate by parts like in the proof of ([70|) as 
many times as we wish. Doing this N times, we obtain 



hi 5, t 



-N 



\O- 3 - N (p)\lp 3 d"\f(p0)f + p 2 d^\f(p0)f}dp 



d6 



+ t 



-N 



\o- 3 - N (p)\[pd^\np0)\ 2 +d^\f(pd)\ 2 }dp 



de 
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oo N-3 



+ t 



—N 



L / £ 10 



2N+l ( 



( P )\d 3 p +l \f(p0)\ 2 dp 



66. 



The top two derivatives of \f\ 2 can be bounded in L 2 , while the remaining lower 
derivatives can be bounded in L°° . Thus, proceeding like in (|70|) . we obtain 



-N 



\\(x) N f\\h + \\(xr- l f\\h + 



N-l . 



N-2 \ 

EiK*>'/iiH 
1=0 / 



<t- N \\(x) N f\\h 



This completes the proof of (|71[) . 



□ 



B.3. Dispersive Decay Estimates. Next, we discuss general L°°-decay estimates 
for the operator H . Since H is of the same form as the linear differential operator 
studied in [19] . we can take advantage of the decay estimate of [HI Theorem 2.2]. 
The proof of this result relies on stationary phase techniques. 

Theorem B.6. Let h = ho, i.e., h(0 — Then, for any k G Z, we have 



(72) 



sup 

igt 3 

sup 

I6E 3 



j>k {Q -ih(£)t+i$-x 

HO 



a fc , 7^, 1 
< f-323(2 fe )-s ) 

1 , 1., 1 



Proof. The first estimate follows directly from [TSJ Theorem 2.2], with phase h and 
amplitude ip k - For the second estimate, we can write 

MO 



HO 



< 



ip k {Oe^ ms+li ' x d^ds 



where we applied the first estimate. The desired bound follows. 

Theorem IB. 61 implies the following dyadic dispersive estimate. 
Corollary B.7. If f e S{R 3 ), t E (0,oo), and < a < \, then 
(73) \\P k H-°e- itH f\\ L ~ <t-i+ a \\P„ k U?f\\ Ll , 

where P^ k = Pk-i + Pk + Pk+i ■ 

Proof. We first consider the case a — 0. Consider the function 



□ 



9k 



9k{x) 



ipk(O e 



-ih(£)t+i£-x 



where h is as in the statement of Theorem lB.6l Then, 



\\P k e- UH f\\ Lr < \\9k*P~kf\\L~ < ||Sik||z-||i > ~fc/IUi. 
The L°°-bound for g k follows from the first estimate of (|72|) . so that 



|P fc e 



-itH 



/Ik? <t-n%{2 k )-^\p^ k f\\ Ll <t-Hp^ k uif\\ L i. 



where in the last step, we applied to the operator U~ x / 2 . 

The case a = 1 is proved similarly, using the second estimate in (|72")l . For the 
remaining < a < 1, the result can be established by interpolation. □ 
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Finally, we piece together the dyadic estimates ([73)) in the various frequency 
bands into a single dispersive estimate. This is stated below, in Proposition lB.8| in 
a form applicable to a general class of Fourier multiplier operators. 

Proposition B.8. Let a',b' € R, and let M be a Fourier multiplier operator, with 
symbol of the form M a i- Also, suppose a € K and I € Z satisfy —3/2 < a < a' , 
I > 0, and b' < I. Then, for any f e S{R 3 ), 

(74) \\e- itH Mf\\ L ~ < m , a >,b>, a ,l i- f ||V'/||^+*- min(l ' 2+a) ||/|U 1 , 

(75) ||e-^M|V|/|| L oo < m , Q))6> ,i t-*||V i+1 /llLi +t- rain( *' 2+a) ||V/|Ui, 
Proof. First, we apply a Littlewood-Paley decomposition to obtain 

\\e~ UH Mf\\ L *o < £ \\P k e- itH Mf\\ L ~ + ]T \\P k e- itH M f\\ L ~ = C + H. 

k<Q k>0 

For H, letting I = b' + e and applying (|68|) and (|69[) . we obtain 

"H < ^2- efe 2( fc ' +£ )' £ ||P fc e- itH /|| L oo < ^2- £k \\P k e- itH V f\\ L ~. 

k>Q k>0 

We can now apply the dyadic dispersive estimate (|73l) to the above. Noting that 
PkU 1 / 2 is bounded on L 1 by (|69p and that 2~ ek is summable yields 

H < t-i 2- £k \\P~kU^ l f\\ L i < HV'/llii- 

k>0 

Next, for the low-frequency terms £, we write a = a' — e, and we apply (|69p to 
M, U 1/2 , and H" " 1 / 2 . This results in the estimate 

C<^2 £k 2 a ' k - £k \\P k e- ltH f\\ L ~ < ^2 £k \\P k e- ltH H^U- 1 * f\\ L ™. 

k<0 k<0 

We can now apply the dyadic dispersive estimate ff73|) . If a < —1/2, then 

c < t- 2 - a 2 e *np~*/iui < *- 2 -»ii/iu*- 

fe<0 

Otherwise, if a > —1/2, then 

£<HE 2 £k \\P^ k H a+ ?f\\ Ll < t-i\\f\\ Ll , 

k<0 

since positive powers of H are bounded on low frequency bands by ([69)) . Combining 
the above bounds for H and C yields ((74)) . 

Finally, for the remaining estimate (fT5"j) . we follow the same steps as above, but 
with / replaced by |V|/, in order to obtain 

He-^AflVl/IU- < *" S E 2- £fe ||P~ fc |V|L^v7|| L1 

k>0 

+ t-iJ2 2ek W p ~k\V\H a+ ?f\\ L1 . 

fe<0 

Using (|6TJ) and ([68]). in addition to <(69]) as before, yields t[75]>. □ 
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B.4. Spatial Decay Estimates. Finally, we obtain some estimates involving spa- 
tial decay of functions. These are used in the proof of Theorem 13.11 in order to 
extract decay rates of subsonic and sonic inertial solutions. Given a multi-index 
a = (at, 012, 013) and any x £ K 3 , we will let x a denote the polynomial x^x^x^ 3 . 

The key component of the decay estimates in the subsonic setting is the following 
Fourier-localized uniform bounds. 

Lemma B.9. Let o,6el, and let M be a Fourier multiplier operator, with symbol 
m in the class A4^- Also, fix a multi-index a — (at, 012,013) and an integer k > 0. 
Then, for any f G <S(R 3 ) and i€l 3 , we have 

(76) \x a P k Mf(x)\ < m , 6 2( 3 / 2 + b > fe ||( a; )H./-|U 2 ,. 

(77) \x a P- k Mf(x)\ < m , a 2(H- Q - 3 ) fe ||( ; r)H/|| L i. 

Proof. First, for (|76j) . by the Fourier inversion formula, we can write 

^■ £ d a [MOm(omm ■ 



\x a P k Mf(x)\ z 

Since i\) k m is in the class MP a and is supported in the band |£| ~ 2 k 



x a P k Mf(x)\< V f \d^f(O\\d 02 bPk(Z)rn(O]\dti 



fi 1 +l3 2 =a 

< 

fl 1 +02=a 

By Holder's inequality, we can control / in L 



\x a P k Mf(x)\<2^ 2+b ^ k \\d p f\\ L > <2 (3/2+b)fc ||(z) H /|| L , 

W<\a\ 

Next, for (l77l). we first write 



\x a P- k Mf(x)\~ / e ix -td a [^ k (0m(li)f(0n 

JM 3 

Since ip- k m is of the form M. h a , 

\x a P^ k Mf(x)\ < V 11^/11^2^1-)* / di 



/3i+/3 2 =a 



< 2(l Q l- a - 3 ) fe ||(x)l Q l/|| L i. □ 

We can construct a decay estimate by summing over the frequency-localized 
estimates of Lemma IB. 91 This following estimate is applied in order to find the 
spatial decay rate for subsonic inertial solutions in Theorem 13. II 



Proposition B.10. Let a£l, and let M be a Fourier multiplier operator, with 
symbol m in the class M.^ 2 ■ If € [0, a + 3), then 

(78) \x\»\Mf{x)\ < m ,^ \\(xr'f\\ L 2, 

for any x 6 M 3 . / G S(R 3 ), and fi' > |>] + 3/2. 
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Proof. We write fj, = fio + fix, where )jlq G Z and G [0, 1). Applying a frequency 
decomposition to the left-hand side of (|78|) and then utilizing Lemma [B. 91 



i \Mf{ X )\ < Y.i^r^iPkMfixw^ixriPkMfix)]} 1 -^ 

k>0 

+Y / [\^r +1 \p-kMf(x)\r[\xr\p^ k Mf(x)\} 1 -^ 

fc>0 

53 2-5 fe || (x) M /|| i2 + ^ 2 ^- a - 3 ) fc || (x> M . 



k>0 k>0 

In particular, we applied (|75|) and (J77J, with \a\ — and |a| = fio + I. Finally, 
since /i — a — 3 < 0, then the above sums converge, and ([78"f follows. □ 

Remark B.ll. One can afeo prove Proposition \B.1(J\ by restricting to the Fourier 
domain |£| > 5, deriving decay estimates independent of 8 in this region, and letting 
5 \ 0. This has the added convenience of avoiding the full dyadic decomposition. 
However, the interpolation between integer powers of decay becomes more compli- 
cated in this case. 

In the sonic setting in Theorem 13. 11 one requires more care in handling the low 
frequencies. This is due to the more singular nature of the operator H~ x near the 
Fourier origin when |u| = 1. The following lemma provides the key component of 
this analysis in the low-frequency region. 

For notational convenience, given a "sonic velocity" ue§ 2 , we define S v to be 
the Fourier multiplier operator given by 



(79) 



(&/)(£) = 



/(0 



In particular, the denominator is especially singular as |£| — > in the u-direction. 

Lemma B.12. Let a,b G R, and let M be a Fourier multiplier operator, with 
symbol m in the class Ai b a . Also, fix v £ § 2 , a multi-index a = (a±, ct2, (J3), and an 
integer k > 0. Then, for any f G <S(R 3 ) and i£l 3 , we have 



(80) 



\x a P- k MS v f(x)\ 



< 



||(a;)l Q l/|| i i2( 3 l"l- Q - 3 ) fc M>0, 
||/|| Ll (l + fc)2-( a + 3 ) fc |a|=0. 



Proof. As in the proof of (|77[) , we begin by writing 



\x a P- k MS v f(x)\ 



^- fc (Qm(Q/(0 

<0 - v ■ i 



and we expand the right-hand side using the Leibniz rule. Derivatives of m, tp-k, 
and / are handled in the same manner as in the proof of ([77]) . For the extra factor, 
we note that if I > 0, then 

VeKit) - * •£)-'] = 9® -(W-vi)- 1 - 1 , 

with g in the class AA°_ 1 . Combining the above, we obtain 



\x a P- k MS v f(x)\ < ^ ||^/lk~2d^l- a ) fe 
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dt 

J 11 f/e\ _ 

< ||(x)l"l/|| L i2(l Q l- a ) fc 



< \\d^f\\ L ~2^ + ^- a ^ k f 



We must now control the integral on the right-hand side. To do this, we convert 
to spherical coordinates (p, 9, 4>) like in the proof of (|2"0j) . with <p € [0, 7r] measuring 
the angle from the w-direction. This yields 



i! 



7/- />2-7T /»7T /■ 2 * 

«S <; / / / P sm ' 



|fi|~a-» ((0 ~ « • 1+H ~ Jo Jo J P ~2-* «P> " cos^+M 

Jo (2-a* + ^)i+N^' 

where in the last step, we used that sin0 < 0, 1— cos</> ~ <fi 2 , and (p)— 1 — p 2 — 2~ 2,£ 
in the domain of integration. This last integral can be evaluated explicitly: 



(2- 2fe + 2 )-H|5 <2 2fc H H>0, 
o (2- 2fc + 2 )i+M "" r ~ 1 log(2- 2 ' £ + 2 )|J < 1 + k \a\=0. 



-dcj) 



< 



Combining the above computations, we obtain (|80|) . □ 

Finally, we state and prove the main estimate pertaining to sonic decay. 

Proposition B.13. Let a £ R, ie< v G § 2 , and ie< M be a Fourier multiplier 
operator, with symbol m in the class .M" 1 . If p £ [0,a/3 + 1), then 

(81) \x\»\MS v f{x)\ <u,n WWfWv, 

for any x £ E 3 , f £ S(R 3 ), and p! > \p] + 3/2. 

Proof. Again, we write p = po + pi, with p £ Z and pi £ [0, 1), so that 
\x\»\MS v f(x)\ <5^[|xr+ 1 |fl k MS„/(a;)|]' l »[|xr|fl fe MS w /(a!)|] 1 -' 41 

fe>0 

+ ^[|xr n + 1 |P-fcAf5 t ,/(x)|]^[| a ;r o |P_ fe A^ t ,/(x)|] 1 -' 11 

fc>0 

= H + L. 

For the high-frequency terms iJ, note that is now in the class A4" 2 , so that, 

like in the proof of (f78|). the estimate (|76[) yields the desired bound for H . For the 
low-frequency terms, we apply 



L < || (x) ^'1 /|| £ i + fc)2[ 3 ' AIn+1 ^ a_3 ^ tlfc 2[ 3 ' t0_a_3 ^ 1_A ' 1 ^ fe 

fe>o 

< ||(a;) M /l|Li + k)2 { ^- a -^ k . 



k>0 

Finally, since 3p — a — 3<0, then the above series converges, and 

L<\\(x)»'f\\ L *. □ 
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